Electronic theory for itinerant in-plane magnetic fluctuations in Na^CoC^ 
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Starting from ab-initio band structure for Na^CoCh, we derive the single-electron energies and the 
effective tight-binding description for the t2 g bands using a projection procedure. We find that due 
to the presence of the next-nearest-neighbor hoppings a local minimum in the electronic dispersion 
close to the F point of the first Brillouin zone forms. Therefore, in addition to a large Fermi surface 
an electron pocket close to the F point emerges at high doping concentrations. The latter yields the 
new scattering channel resulting in a peak structure of the itinerant magnetic susceptibility at small 
momenta. This indicates itinerant in-plane ferromagnetic state above certain critical concentration 
Xm, m a greement with neutron scattering data. Below x m the magnetic susceptibility shows a 
tendency towards the antiferromagnetic fluctuations. We estimate the value of 0.56 < x m < 0.68 
within the rigid band model and within the Hubbard model with infinite on-site Coulomb repulsion 
consistent with the experimental phase diagram. 

PACS numbers: 31.15.Ar; 74.70.-b; 71.10.-w; 75.40.Cx 



1. The diverse physical properties of the cobaltate 
Na^CoC^ attracted much attention after the discovery 
of the unconventional superconductivity in its hydrated 
counterpart, Na^CoC^ • 2/H2O The phase diagram of 
this compound, with varying electron doping x and water 
intercalation y, is rich and complicated; in addition to su- 
perconductivity, it exhibits magnetic and charge orders, 
and some other structural transitions 0, S, 0, @| ■ Par- 
ent compound, Na T Co02, is a quasi-two-dimensional sys- 
tem with Co in C0O2 layers forming a triangular lattice 
where Co-Co in-plane distance three times smaller than 
the inter-plane one. Na ion resides between the C0O2 lay- 
ers and gives additional x electrons to the layer, lowering 
Co valence from Co 4+ (3d 5 configuration) to Co 3+ (3d 6 
configuration) upon x changing from for virtual com- 
pound C0O2 to 1 for NaCo02- The hole in the d-orbital 
occupies one of the t^g levels, which are lower than e g 
levels by about 2 eV [6| . The degeneracy of the t^ g levels 
is partially lifted by the trigonal distortion which splits 
it into the higher a\ g singlet and the lower two e' g states. 

First principle LDA (local density approximation) and 
LDA+U band structure calculations predict Na K Co02 
has a large Fermi surface (FS) having mainly a\ g char- 
acter and centered around the T = (0, 0, 0) point and 
also six hole pockets of mostly e' g character near the 
K = (0,-^,0) points of the hexagonal Brillouin zone 
for a wide range of x [1, Q. At the same time, recent 
Angle-Resolved Photo-Emission Spectroscopy (ARPES) 
experiments 

HIE , [f0 | reveal doping dependent Fermi 
surface evolution for a wide range of Na concentrations 
(0.3 < x < 0.8) with no sign of the hole pockets. The ob- 
served Fermi surface is centered around the T point and 
have mostly a\ g character. Furthermore, measured dis- 
persion of the top of the valence band is twice as narrower 



as compared to the LDA calculated bands. 

Concerning magnetic properties the local spin den- 
sity approach (LSDA) predicts Na x Co02 to have a weak 
intra-plane itinerant ferromagnetic (FM) state for almost 
all Na concentrations, 0.3 < x < 0.7 [12j. On the con- 
trary, neutron scattering finds the A-type antiferromag- 
netic order implying the ferromagnetic order within Co- 
layer only for 0.75 < x < 0.9 with ordering temperature 
T m « 22 K, and with inter-plane J c and intra-plane J a b 
exchange constants to be 12 meV and -6 meV, respec- 
tively [iiiiiiii- 

In this letter we derive an effective low-energy model 
describing the bands crossing the Fermi level on the basis 
of the LDA band structure calculations. Due to the FS 
topology, inferred from LDA band structure, the mag- 
netic susceptibility xo (cu to = 0) reveals two different 
regimes for different dopings: for x < 0.58 it shows pro- 
nounced peaks at antiferromagnetic (AFM) wave vector 

Qafm = {(t> 73) > (t> )} resultm S in the tendency 
towards the in-plane 120° AFM order, while for x > 0.58 
the peaks at small momenta near the Qfm = (0, 0) form, 
clearly demonstrating the tendency of the system towards 
the itinerant in-plane FM ordered state. We find that 
the formation of the electron pocket around the T point 
is crucial for the in-plane FM ordering at high doping 
concentrations . 

2. The band structure of Na .33CoO 2 (Fig. [TJ was 
obtained within the LDA [l6| in the framework of TB- 
LMTO-ASA (Tight Binding approach to the LMTO 
using Atomic Sphere Approximation) [ijj computation 
scheme. The crystal structure parameters were taken 



from [18[ . The bands crossing the Fermi level have mostly 
dig character, consistent with previous LDA findings [fj]. 
Note, the small FS pockets near the K point with e'- 
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FIG. 1: Calculated near-Fermi level LDA band structure and 
partial density of states (PDOS) for Nao.33Co02. The contri- 
bution of Co-aig states is denoted by the vertical broadening 
of the bands with thickness proportional to the weight of the 
contribution. The crosses indicate the dispersion of the bands 
obtained by projection on the ti g orbitals. 



symmetry present at x = 0.33 disappear for higher dop- 
ings because the corresponding bands sink below the 
Fermi level. 

To construct the effective Hamiltonian and to derive 
the effective Co-Co hopping integrals t°f for the Un- 



fg 

manifold we apply the projection procedure [191 [2CJ. 
Here, (a/3) denotes a pair of orbitals, a\ g , e' gl or e' g2 - 
The indices / and g correspond to the Co-sites on the 
triangular lattice. The obtained hoppings are given in 
the Table [IJ and the obtained single-electron energies 
e a are equal to (in eV, relative to £ ai '): e aig = 0, 
e<i = e e '^ = -0.053. 

An agreement between the bands obtained using pro- 
jection procedure and the LDA bands shown in Fig. [TJ 
that confirms the Co-i2g nature of the near-Fermi level 
bands @, HU. For simplicity we have enumerated site 



pairs, t 



a/3 
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t%P, with n = 0,1,2,... (see Fig. Hi and 
the correspondence between in-plane vectors and index 
n in the Table []}. Due to C3 symmetry of the cobaltate 



lattice, the following equalities arc present: \t 
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a\ g — ► aig hoppings, which, however, does not hold for 
e g i 2 orbitals. Thus, since the hybridization between the 
aig and the e' g bands is not small, a simplified descrip- 
tion of the bands crossing the Fermi level in terms of the 
a,\ g band only (neglec ting e' g band and corresponding hy- 
bridizations, see e.g. [221 ]) may lead to a incorrect result 
due to its higher symmetry. In the following we neglect 
the inter-layer splitting present for k z = plane because 
of its subtle effect on the topology of the FS [f| . 

Then the free electron Hamiltonian for Co02-plane in 
a hole representation is given by: 



EEcc^,. (i) 



k,a,cr 



k,fr a,/3 



where a\ aa (d\^ aa ) is the annihilation (creation) operator 




FIG. 2: (a) Schematic crystal structure of the Co-layer in 
Naj;Co02 with hopping notations within the first three co- 
ordination spheres (C.S.). (b) LDA-calculated Fermi surface 
with cylindrical part having mostly a\ g character and six hole 
pockets having mostly e' g character. k x and k y coordinates of 
the symmetry points are given in units of 2n/a with a being 
the in-plane lattice constant. 



hopping matrix element 



for hole with momentum k, spin a and orbital index a, 
= rfit„„rfka<Tj and i£ i s the Fourier transform of the 
Introducing matrix notations, 

/a/3 ~ "k a ^ ynJafj = l n P > the no PP in S s matrix 
elements in the momentum representation are given by: 



tZ and (in 



£k = 2ti cos &2 + 2i2 cos &3 + 2^3 cos ki + 

+ 2t4eos(/ci + £3) + 2t5 cos(/c2 + ki) + 2ig cos(fci — fe) + 
+ 2t 7 cos 2fc 2 + 2t 8 cos 2A: 3 + 2t 9 cos 2fci , (2) 



where k\ = ^k 



Within this rigid band approximation the doping- 
dependent evolution of the electronic dispersion, the den- 
sity of states (DOS) and the FS is shown in Fig. [3] We 
notice, that already at x = 0.48 the FS e' g hole pockets 
are absent. Most importantly we find another interest- 
ing feature. Namely, the local minimum of the band dis- 
persion around the T point yields the appearance of the 
second FS contour centered around this point. This elec- 
tron FS pocket becomes larger upon increasing doping x. 
As was shown in Ref. [23| for the Hubbard model on a 
triangular lattice, the main reason for the local minimum 
around the T point is the presence of the next-nearest- 
neighbor hoppings, which also enter in our calculations. 
Although this minimum is not yet directly observed in 
ARPES experiments, the presence of the associated sec- 
ond FS contour would reduce the FS volume and resolve 
an issue why the volume of the FS observed in ARPES is 
larger than it is expected from Luttinger's theorem 24 1. 
Furthermore, an emergence of this pocket would influ- 
ence the Hall conductivity at high doping concentrations 
which is interesting to check experimentally. 

3. To analyze a possibility of the itinerant magnetism 
we calculate the magnetic susceptibility xo(q,u = 0) 
based on the Hamiltonian Hq. The doping-dependent 
evolution of the peaks in Re%o (q, 0) is shown in Fig. [4] 
At x = 0.45 the e' g bands are below the Fermi level, and 
the FS has the form of the rounded hexagon. It results in 
a number of nesting wave vectors around the antifcrro- 
magnetic wave vector Qafai- The corresponding broad 
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TABLE I: The in-plane hopping integrals t„ for different in-plane vectors n = (/, g) for Na^CoO^, x = 0.33. 
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FIG. 3: Calculated band structure and Fermi surface topology 
for Na^CoO-2 for x = 0.47, 0.68, 0.73. The dashed (light gray) 
and solid (black) curves represent the rigid-band approxima- 
tion and the Hubbard-I solution, respectively. The horizontal 
line denotes self-consistently calculated chemical potential [i. 



peaks in the Rexo(QiO) appear around Qafm, indicat- 
ing the tendency of the electronic system towards the 
120° AFM ordered state [25|. Upon increasing doping, 
the Fermi level crosses the local minimum at the V point, 
resulting in the second almost circle FS contour. As soon 
as this change of the FS topology takes place, the scat- 
tering at the momentum Qafm is quickly suppressed, 
and vanishes already at x m ~ 0.56. Most importantly, a 
new scattering vector, Qi, appears. This wave vector is 
small and yields peaks in the magnetic susceptibility at 
small momenta, indicating the tendency of the magnetic 
system to shift towards itinerant FM order. For larger 
x the inner FS contour increases leading to a further de- 
crease of the Qi. In the case of i k 0.88 the FS topol- 
ogy changes again, resulting in a six distant FS contours, 
yielding even smaller length of Qi. The obtained scat- 
tering at small momenta in the magnetic susceptibility 
for x > x m is qualitatively consistent with the scattering 
around Qfm = (0, 0), observed in the neutron scattering 
experiments 13, 1J, 15j |. 

4. Since obtained magnetic susceptibility depends 




FIG. 4: The grayscale plot of the real part of the magnetic 
susceptibility Rexo(k-,LJ — 0) as a function of the momentum 
in units of 27r/a (left), and the Fermi surface for corresponding 
doping x (right). The arrows indicate the scattering wave 
vectors Qi as described in the text. 



mostly on the topology of the FS one expects that the 
behavior shown in Fig. |4] will be valid even if one includes 
an interaction term Hi nt into account, at least in the case 
if it is the on-site Hubbard interaction U. The only dif- 
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ference would be a shift of the critical concentrations x m , 
at which the FS topology changes and tendency to the 
AFM order changes towards the FM ordered state. To 
check this, we have taken the strong electron correlations 
into account by adding the on-site Coulomb interaction 
terms to the Ho, similar to Refs. [26|, 27 1. The effec- 
tive on-site Hubbard repulsion U e fj w 4 eV on Co sites 
is much larger than the bare bandwidth W « 1.2 eV, 
and, thus, it is possible to project doubly occupied states 
out and formulate an effective model equivalent to the 
Hubbard model with infinite value of U. 

In the atomic limit local low-energy states on the Co 
sites are the vacuum state |0) and the single- hole states 
aer), |ei<r), |e2<x). The simplest way to describe the 
quasiparticle excitations between these states is to use 



the projective Hubbard AT-operators [28|: X 
\p) (q\, where index m <-> (p,q) enumerates quasipar- 
ticles. There is a simple correspondence between X- 
operators and single-electron creation-annihilation oper- 
ators: dfaa = laai^XJ 1 , where 7 acr (m) determines 

the partial weight of a quasiparticle m with spin a and 
orbital index a. In these notations the Hamiltonian of 
the Hubbard model in the limit U — > oo has the form: 



H 



f,p 



E E t% m Xfxf. (3) 

f=£g m,m' 



In the so-called Hubbard-I approximation within the gen- 
eralized Dyson equation for the V-operators 2^, 3^, 31 1 
the quasiparticle bands formed by the a\ g — > a\ g hop- 
pings will be renormalized by the (1 + x)/2 factor, while 
the quasiparticle bands formed by the e' g hoppings will 
be renormalized by x. 

In Fig. [3] the quasiparticle spectrum is shown. One 
finds within Hubbard-I approximations the bands be- 
come narrower with lowering x due to doping dependence 
of the quasiparticle's spectral weight. Most importantly, 
the doping evolution of the FS is qualitatively the same 
as in the rigid-band picture. Thus, a bandwidth reduc- 
tion and a spectral weight renormalization do not change 
the topology of the FS. Therefore, the presence of the 
strong electronic correlations do not change qualitatively 
our results for the bare susceptibility. Quantitatively, the 
critical concentration x m shifts towards higher values and 
within Hubbard-I it becomes x m w 0.68. 

5. To conclude, we have shown that in the model with 
ab-initio calculated parameters the magnetic susceptibil- 
ity is doping dependent. At the critical doping concen- 
tration, x m , the electron pocket on the FS in the center 
of the Brillouin zone well develops. For x < x m , the 
system shows tendency towards the 120° AFM ordered 
state, while for x > x m the peak in the magnetic suscepti- 
bility is at small wave vectors indicating strong tendency 
towards the itinerant FS state. Within the tight-binding 
model x m is estimated to be around 0.56. Analyzing the 
influence of the strong Coulomb repulsion and the corre- 



sponding reduction of the bandwidth and the quasipar- 
ticle spectral weight in the strong-coupling Hubbard-I 
approximation, we show that the critical concentration 
changes to the x m « 0.68. At the same time, the under- 
lying physics of the formation of the itinerant FM state 
remains the same. 
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